NASA Technical Paper 3466 



Application of Two Procedures for Dual-Point 
Design of Transonic Airfoils 

Raymond E. Mineck, Richard L. Campbell, and Dennis O. Allison 
Langley Research Center • Hampton, Virginia 


National Aeronautics and Space Administration 
Langley Research Center • Hampton, Virginia 23681-0001 


September 1994 



This publication is available from the following sources: 


NASA Center for AeroSpace Information 
800 Elkridge Landing Road 
Linthicum Heights, MD 21090-2934 
(301) 621-0390 


National Technical Information Service (NTIS j 
5285 Port Royal Road 
Springfield, VA 22161-2171 
(703) 487-4650 



Abstract 


Two dual-point design procedures were developed to reduce the objective 
function of a baseline airfoil at two design points. The first procedure to 
develop a redesigned airfoil used a weighted average of the shapes of two 
intermediate airfoils redesigned at each of the two design points. The second 
procedure used a weighted average of two pressure distributions obtained from 
an intermediate airfoil redesigned at each of the two design points. Each 
procedure was used to design a new airfoil with reduced wave drag at the 
cruise condition without increasing the wave drag or pitching moment at 
the climb condition. Two cycles of the airfoil shape- averaging procedure 
successfully designed a new airfoil that reduced the objective function and 
satisfied the constraints. One cycle of the target (desired) pressure- averaging 
procedure was used to design two new airfoils that reduced the objective 
function and came close to satisfying the constraints. 


Introduction 

Transport aircraft are typically designed to cruise 
efficiently at a single flight condition, the cruise de- 
sign point, while meeting other design constraints. 
As the result of a myriad of factors such as air traf- 
fic control operations, flights behind schedule, and 
the change in weight as fuel is burned, or other fac- 
tors, aircraft are often operated at less efficient, off- 
design flight conditions. The penalty for operating 
at off-design flight conditions can be significant, es- 
pecially for wings with supercritical airfoils. Some 
aircraft, such as the high-speed civil transport, may 
have two cruise design flight conditions: one for op- 
eration over land and the other for operation over 
water. Aerodynamic performance may be improved 
if a dual-point design procedure can be developed 
that considers two design points. Because the steps 
of a dual- point design procedure for a wing should 
be similar to those for an airfoil, dual-point design 
procedures can first be developed for the simpler 
two-dimensional (2-D) problem and then be extended 
to the three-dimensional (3-D) problem. Two dual- 
point design procedures arc proposed herein to re- 
design an airfoil for reduced wave drag at two design 
points. 

Airfoil design codes have been developed that de- 
fine the airfoil associated with a specified pressure 
distribution at a single design point. The codes fall 
into two general categories: inverse methods and 

predictor/corrector methods. Inverse methods, such 
as the method described in reference 1, determine 
the airfoil shape by directly using the velocity dis- 
tribution on the airfoil surface as a boundary condi- 
tion. The velocity distribution is derived from the 
specified surface pressure distribution. Predictor/ 
corrector methods, such as the method described in 


reference 2, determine the airfoil shape in an iterative 
process. The predictor module computes the flow 
field and pressure distribution associated with the 
current airfoil shape. The corrector module adjusts 
the airfoil shape to bring the computed and specified 
pressure distributions into agreement. To ensure that 
the airfoil shape is practical, constraints on the shape 
are often included in the corrector module, which can 
be a numerical optimization technique such as the 
technique described in reference 3. Numerical opti- 
mization techniques search for the minimum value 
of an objective function subject to constraints. Use 
of numerical optimization entails a significant in- 
crease in the computational time because of the large 
number of cases that must be evaluated. 

Dual-point design capability can be added to 
predictor/corrector methods that use numerical op- 
timization by adding constraints on the appropriate 
parameters at the second design point or by including 
the second design point in the objective function to 
be minimized. Application of numerical optimization 
for the dual-point design of an airfoil was demon- 
strated in reference 4. The goal was to reduce the 
drag at the maximum operating Mach number with- 
out increasing the drag at the maximum range Mach 
number. The airfoil thickness distribution was con- 
strained such that the original profile had to be con- 
tained within the redesigned profile. A simple objec- 
tive function to minimize the drag at the maximum 
operating Mach number was defined. The drag at 
the maximum range Mach number was constrained 
to be less than that of the baseline airfoil. The re- 
sults indicated that the goal was achieved. However, 
satisfying the constraints at the second design point 
reduced the benefits obtained at the first design point 
without considering the constraint. 



Numerical optimization was used in reference 5 to 
reduce the drag on an airfoil at two flow conditions. 
In this case, the objective function was the weighted 
sum of the drag coefficients at each flow condition. 
No constraints were imposed on the pitching moment 
of the redesigned airfoil. Note that the optimum 
shape is for the set of shape functions used to re- 
design the airfoil shape and that the optimum shape 
may not represent the global optimum. Dual-point 
designs using numerical optimization require a large 
number of geometries to be evaluated. The num- 
ber of airfoil geometries that were evaluated in the 
dual-point design process for the cases studied ranged 
from 61 to 105. The associated computer resources 
required to evaluate these geometries is large com- 
pared with the resources required for a single-point 
design without numerical optimization. 

Alternative procedures, which require fewer com- 
puter resources than current numerical optimization 
techniques, are desired to redesign an airfoil at more 
than one design point. A simple multipoint design 
procedure was used in reference 6 to improve the 
performance of a cascade at off-design conditions. 
The procedure analyzed the baseline cascade airfoil 
at angles of attack from —3° to 3° from the de- 
sign operating condition to obtain the pressure dis- 
tributions at these two off-design conditions. The 
two pressure distributions were modified to obtain 
improved boundary- layer characteristics, and these 
modified pressure distributions were then used to de- 
sign two intermediate airfoils, one at each off-design 
point. The pressure distribution at the design op- 
erating condition, which was computed for each of 
the two intermediate airfoils, was then used to de- 
velop a new target (desired) pressure distribution at 
the design point. This new target pressure distribu- 
tion was used to design a new airfoil at the design 
operating condition. This procedure required seven 
analysis cases and three single-point design cases, sig- 
nificantly fewer than the number of cases processed 
for the numerical optimization study discussed in ref- 
erence 5. Results indicated that the new airfoil had 
better performance than the original airfoil over an 
angle-of-attack range from —2° to 2° from the design 
operating condition. Although the new airfoil may 
not be the global optimum, the significant reduction 
in computer resources makes this type of procedure 
attractive. 

Two dual-point design procedures were proposed 
in reference 7 that did not require a large number 
of cases to be processed to define a new airfoil. The 
two design points arc assumed to be sufficiently close 
to each other that the aerodynamic characteristics 
at the two design points are similar to each other. 


Changes in the aerodynamic characteristics are as- 
sumed to vary linearly with changes in the airfoil 
shape. This linear assumption is the underlying basis 
for the two dual-point design procedures. One proce- 
dure used the simple average of the shapes of two in- 
termediate airfoils, with one designed at each design 
point, to define a new airfoil. The other procedure 
used the simple average of pressure distributions of 
two intermediate airfoils to design a new airfoil. Each 
procedure used the same predictor /corrector airfoil 
design code. The corrector module included options 
to constrain the values of certain aerodynamic and 
geometric characteristics. These procedures required 
the user to define the “optimum” pressure distribu- 
tion to be used in the dual-point design process at 
each design point. 

Under a cooperative agreement between NASA 
and the Cessna Aircraft Company, the dual-point de- 
sign procedures from reference 7 were used to design 
a new airfoil, and wind tunnel tests were conducted 
to verify the results. With user assistance, the two 
dual- point design procedures were used sequentially 
to define a redesigned airfoil that had less wave drag 
at each design point subject to constraints on the 
pitching moment and the airfoil geometric character- 
istics. Models of the baseline and the redesigned air- 
foils were tested in the Langley 0.3-Meter Transonic 
Cryogenic Tunnel (0.3-m TCT) described in refer- 
ence 8. Results showed total drag reductions slightly 
larger than the predicted wave-drag reductions. De- 
tails of the tests of the baseline airfoil are presented 
in reference 9. 

Information about the dual-point design proce- 
dures, the development of the redesigned airfoil, and 
the experimental verification of the reduced drag of 
the redesigned airfoil is included in the appendix 
of this report. The results indicate that averag- 
ing is a viable technique for use in a dual- point de- 
sign procedure. The good agreement between the 
computed results and the experimental results in- 
dicates that the Euler flow solver with the interac- 
tively coupled boundary layer is a suitable mathe- 
matical model for the flow fields encountered in this 
design process. However, the dual- point design pro- 
cedures of reference 7 had two major drawbacks: the 
user must define “optimum” pressure distributions, 
and the simple averaging technique did not necessar- 
ily produce the best design. This report describes 
revised dual-point procedures that eliminate these 
drawbacks. 

The purpose of this report is to describe two re- 
vised dual-point design procedures and to demon- 
strate their application to redesign an airfoil to re- 
duce the value of an objective function specified by 
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the user. These revised dual-point design procedures, 
which require significantly fewer cases to be pro- 
cessed than procedures using numerical optimization, 
should require fewer computer resources. This report 
presents a description of two revised dual-point de- 
sign procedures and the associated technique for con- 
strained airfoil design. Application of each modified 
procedure to a climb condition and to a cruise con- 
dition is presented as a practical demonstration for 
different objective functions. 

Symbols 

The results are presented in coefficient form with 
the moment reference center at the quarter-chord. 
All experimental measurements and calculations 
were made in U. S. customary units. 

C p local static pressure coefficient, P q^°° 

C* sonic pressure coefficient 

c chord, ft 

Crf section drag coefficient, 

A Cfj change in section drag coefficient 

, . , re ■ i. Wave drag 

°d w section wave-drag coefficient, q^c — 

C[ section lift coefficient, 

Cm section pitching-moment coefficient, 

Pitching moment 

<hc <- 1 

F objective function 

G constraint function 

Moo free- stream Mach number 

p local static pressure, lbf/ft 2 

p 0 c free-stream static pressure, lbf/ft 2 

free-stream dynamic pressure, 
ibf/ft 2 

R c Reynolds number based on chord length, 

Po cVqc c 
Poe 

ri e leadingTedge radius, ft 

t local thickness, ft 

VJx, free-stream velocity, ft /sec 

W weighting factor 

x chordwise distance, positive measured aft 

from leading edge, ft 


y vertical distance, positive measured up 

from airfoil reference line, ft 

a angle of attack, positive leading edge 

up, deg 

/ioc free-stream viscosity, slugs/ft-sec 

Poc free-stream density, slugs/ft 3 

Subscripts: 

/ data at foot of shock 

max maximum 

s shock 

1 design point 1 

2 design point 2 
Abbreviations: 

DISC Direct Iterative Surface Curvature 
TCT Transonic Cryogenic Tunnel 

Airfoil Constrained Design Technique 

The constrained design technique of reference 10 
was selected for use in this study. This technique 
designs a new airfoil that produces the desired pres- 
sure distribution at the specified flow condition (de- 
sign point) subject to user-defined constraints. In 
this study, the new airfoil was subject to the same 
geometric and aerodynamic constraints that were ap- 
plied to the design of the baseline airfoil. The tech- 
nique is divided into three primary parts: a con- 

strained design module to adjust the target (desired) 
pressure distribution to satisfy the constraints, an 
airfoil design module to modify the airfoil shape to 
achieve the target pressure distribution, and an aero- 
dynamic analysis module (flow solver) to compute 
the pressure distribution associated with the current 
airfoil shape. These three parts are used in an iter- 
ative process to determine the new airfoil shape as 
shown in figure 1. 

The initial airfoil shape and the initial target pres- 
sure distribution are input into the computer pro- 
gram to define the current airfoil shape and target 
pressure distribution. The aerodynamic analysis 
module then computes the pressure distribution as- 
sociated with the current airfoil shape. The airfoil 
design module adjusts the airfoil shape after each 
iteration to bring the current computed pressure dis- 
tribution into agreement with the current target pres- 
sure distribution. After a specified number of it- 
erations, the constrained design module compares 
the geometric and aerodynamic characteristics of the 


3 



current airfoil shape with the constraints. If any con- 
straints are violated, the constrained design module 
adjusts the target pressure distribution accordingly. 
The process is repeated for a user-specified number 
of design iterations. Efficiency is improved over se- 
quential techniques because the airfoil shape and tar- 
get pressure distribution converge concurrently. A 
brief description of the three primary parts of the 
constrained design technique follows. 

The constrained design module modifies the de- 
tailed target pressure distribution to eliminate viola- 
tions of the user-specified constraints. Details of the 
module are presented in reference 10. The detailed 
target pressure distribution consists of the target (de- 
sired) pressure coefficient at each specified chord wise 
station. Seven target control points along the chord 
of each surface of the airfoil are determined to define 
the major characteristics of the detailed target pres- 
sure distribution. Seven was the smallest number of 
control points required to define reasonably well the 
acceleration, rooftop, shock, and deceleration regions 
of the chordwise pressure distribution. A sample of 
a detailed target pressure distribution and the asso- 
ciated target control points is presented in figure 2. 
These target control points are adjusted by the con- 
strained design module to satisfy the aerodynamic 
and geometric constraints. The detailed target pres- 
sure distribution is adjusted to match the adjusted 
target control points. The initial target pressure dis- 
tribution can be identical to the pressure distribution 
obtained from the aerodynamic analysis module for 
the baseline airfoil shape at the specified design point 
(Mach number, Reynolds number, and section lift co- 
efficient). The user may optionally specify the initial 
target pressure distribution to achieve specific goals 
such as reduced adverse pressure gradient or reduced 
Mach number ahead of the shock. 

The constrained design module monitors both 
the aerodynamic and the geometric constraints. As 
used in this project, the aerodynamic constraints are 
the lift coefficient, wave-drag coefficient, pitching- 
moment coefficient, and the pressure gradient on 
both surfaces over the forward portion of the air- 
foil. The geometric constraints are the leading-edge 
radius, maximum thickness, and local thickness at 
a specified chordwise location. The constrained de- 
sign module computes the required aerodynamic and 
geometric characteristics from the current pressure 
distribution and from the current airfoil shape. If 
any constraints are violated, the target control points 
are adjusted based on empirically derived expressions 
described in reference 10. In certain instances, the 
expressions may indicate that a fully constrained tar- 
get is not possible. In these cases, one or more aero- 


dynamic constraints are relaxed to allow the remain- 
ing constraints to be satisfied. The pitching-moment 
coefficient is relaxed first and the wave-drag coeffi- 
cient is relaxed second. The adjusted target control 
points are used to adjust the detailed target pressure 
distribution. 

The airfoil design module utilizes the Direct 
Iterative Surface Curvature (DISC) method of 
Campbell and Smith (ref. 2) which designs a new 
airfoil by iterating the airfoil shape to obtain the 
target pressure distribution. The airfoil design mod- 
ule compares the computed current pressure distri- 
bution associated with the current shape with the 
detailed target pressure distribution. The difference 
between the computed current and target pressure 
distributions at each chordwise location on the air- 
foil is used to determine the change in the local airfoil 
surface curvature. The ordinates at all chordwise lo- 
cations downstream of the curvature change are lin- 
early sheared to obtain the required change in local 
curvature, as shown in figure 3(a). If the resulting 
airfoil does not close at the trailing edge, the points 
on each surface are linearly sheared to form a conven- 
tional, closed trailing edge, as shown in figure 3(b). 
The nose camber line and the airfoil upper and lower 
surfaces are smoothed. The computational grid is 
then automatically adjusted to compensate for the 
change in the airfoil shape. 

Since this application could involve flow fields 
with strong shocks, the Euler equations were se- 
lected to model the flow field. The Euler equation 
solver GAUSS2, which was developed by Hartwich 
(ref. 11), was used for the aerodynamic analysis mod- 
ule. code solves the Euler equations, cast in 

a nonconservative form, on a structured grid using 
an implicit upwind-differencing procedure. The use 
of a shock-fitting approach produces sharp shocks, 
even on relatively coarse grids. In regions away from 
the shock, a second-order, split-cocfficient-matrix up- 
winding method is used. Convergence to steady state 
is accelerated by a diagonalized approximate fac- 
torization technique. The airfoil lift and pitching- 
moment coefficients are computed by integrating the 
chordwise pressure coefficient distribution. Viscous 
effects have been incorporated into the code by using 
a boundary-layer displacement-thickness approach. 
For this study, the boundary-layer characteristics are 
computed in a separate? subroutine by using a modi- 
fied version of the approach by Stratford and Beavers 
(ref. 12), and the viscous drag coefficient is estimated 
by using the technique by Squire and Young (ref. 13). 
Wave drag is computed from the far-field entropy, 
based on the approach of Oswatitsch (ref. 14). 
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Description of Revised Dual-Point 
Design Procedures 

The two dual- point design procedures of refer- 
ence 7 (also included in the appendix of this re- 
port) were revised for use in this study. The revi- 
sions eliminated the need for the user to specify an 
initial improved target pressure distribution and de- 
fined a weighted averaging method based on reduc- 
ing the objective function. The goal of each proce- 
dure was to define a redesigned airfoil that, subject 
to constraints, had a smaller objective function than 
the baseline airfoil. The objective function could in- 
volve aerodynamic and geometric characteristics of 
the airfoil at more than one design point. Possi- 
ble constraints include the airfoil leading-edge radius, 
maximum airfoil thickness, thickness at a chordwise 
station, pitching moment, and wave drag at either of 
the design points. 

Initial Processing for Each Design Cycle 

The first dual-point design procedure uses the 
weighted average of the ordinates of two inter- 
mediate airfoils to define a new airfoil shape, and 
this procedure will be referred to as “airfoil shape 
averaging.” The second dual-point design procedure 
uses the weighted average of two pressure distribu- 
tions to design a new airfoil, and this procedure will 
be referred to as “target pressure averaging.” Each 
procedure begins with the same initial processing, as 
shown in figure 4. The baseline airfoil A is analyzed 
at both design points. The desired aerodynamic and 
geometric characteristics are specified at each design 
point for use by the airfoil constrained design tech- 
nique, that is, the constrained DISC method. This 
method modifies the pressure distribution from the 
analysis of the baseline airfoil to achieve the desired 
aerodynamic and geometric characteristics so that 
the user is not required to define an improved ini- 
tial target pressure distribution, as was the case in 
reference 7. Two intermediate airfoils, designated B 
and C, are designed, one at each design point, by 
using the airfoil constrained design technique. The 
aerodynamic characteristics of each intermediate air- 
foil, developed at one design point, are analyzed at 
the other design point. By using these results, the 
objective function ( F ) is evaluated for the baseline 
airfoil and for the two intermediate airfoils. 

The objective function is plotted against the 
weighting factor of one of the airfoils (VP), as shown 
in figure 5. For this example, the objective function 
for airfoil B is plotted at 1.0 (all of airfoil B and none 
of airfoil C), and the objective function for airfoil C 
is plotted at 0 (none of airfoil B and all of airfoil C). 
The value of the objective function is assumed to 


vary linearly with the fraction of the airfoil. Con- 
straint functions ( G ) are also plotted as a fraction of 
the same airfoil. If none of the constraints are vio- 
lated, the intermediate airfoil, either B or C, having 
the lower objective function is the new airfoil. If one 
or more of the constraints are violated, the weighting 
factor of airfoil B (fig) that produces the minimum 
objective function without violating the constraints is 
selected. This value of the weighting factor is used for 
the procedures for airfoil shape averaging and target 
pressure averaging. 

Use of the weighting factor results in selecting 
the smallest value of the objective function for which 
no constraint violations occur, thus eliminating the 
drawback of the simple average used in reference 7. 
In the initial processing, the airfoil constrained de- 
sign technique developed an intermediate airfoil at 
each design point that had improved characteris- 
tics relative to the baseline airfoil at that particu- 
lar design point. Changes to the airfoil shape or 
to the pressure distribution are associated with the 
improved characteristics, and the weighted average 
serves as a filter. Similar changes at a chordwise lo- 
cation for both intermediate airfoils will pass through 
the weighted average and promote those changes in 
the new airfoil. Opposing changes will cancel part 
of each other and reduce those changes in the new 
airfoil. Thus, the new airfoil shape or the new tar- 
get pressure distribution from the weighted average 
should contain some of the improvements designed 
into the intermediate airfoils by the constrained 
airfoil design technique. 

Airfoil Shape Averaging 

The airfoil shape- aver aging procedure uses a 
weighted average of the shapes of the two interme- 
diate airfoils to define a new airfoil (D), as shown in 
figure 4. The procedure begins with the initial pro- 
cessing described previously, shown in the left part 
of figure 4, which defines two intermediate airfoils. 
The aerodynamic and geometric characteristics are 
used to define the values of the objective function 
and the constraints. These are then used to obtain 
the weighting factor (VVpj). The new r (redesigned) air- 
foil D is the weighted average of the shapes of tw r o 
intermediate airfoils, B and C. At each chordwise 
location, the ordinate of the redesigned airfoil D is 
defined by 

-(D)= Wb -(B) + (1 - W 7 b) -(C) 

c c c 

By using the weighted average, some improvements 
from the constrained DISC method are passed to the 
averaged airfoil. The aerodynamic characteristics of 
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the averaged airfoil D are analyzed at both design 
points. The results are used to evaluate the objective 
function and to check for any constraint violations. If 
the objective function for the averaged airfoil is less 
than the objective function for the baseline airfoil 
and no violated constraints occur, the dual-point de- 
sign cycle is complete. To find an airfoil shape closer 
to the optimum, the new airfoil D becomes the base- 
line airfoil A. and then the airfoil shape- averaging 
procedure is repeated by beginning with the initial 
processing. If the objective function for the new (av- 
eraged) airfoil is greater than the objective function 
for the baseline airfoil, different intermediate airfoils 
are needed. Less stringent aerodynamic and/or geo- 
metric' constraints should be specified at each design 
point for use by the airfoil constrained design tech- 
nique. Experiences with the procedures suggest that 
reducing the desired reduction in the wave drag by 
about half will often alleviate problems in obtaining 
satisfactory intermediate airfoils. The initial process- 
ing is then repeated to define two different intermedi- 
ate airfoils. The airfoil shape-averaging procedure is 
repeated by using the different intermediate airfoils. 

Airfoil Target Pressure Averaging 

The target pressure- averaging procedure uses a 
weighted average of the pressure distributions of the 
two intermediate airfoils at one of the design points to 
define a new target pressure distribution, as shown in 
figure (i. The procedure begins with the initial pro- 
cessing described previously (see left part of fig. 6) 
that defines two intermediate airfoils. The aero- 
dynamic' and geometric characteristics are used to 
define the value's of the objective function and the 
constraints. These are then used to obtain the 
weighting factor (Up). The desired pressure distri- 
bution for t he new (redesigned) airfoil, either E or F, 
is the' weighted average of the chord wise pressure dis- 
tributions of the two intermediate airfoils, B and C. 
At each chordwise location, the desired target pres- 
sure 4 distribution of the redesigned airfoil is defined 
for the first design point by 

C llA (E)=W li C pA (B)+(l - Hr) Cp.i(C) 

and for the second design point by 

C /j .2(F)= WbC,.j(B) + (1 - Wb) <7*2(0 

By using a weighted average of these tw r o interme- 
diate target pressure distributions, effects from the 
design of both intermediate airfoils are included in 
the new 7 target pressure distribution. 

The new target pressure distribution is the 
weighted average of the intermediate pressure distri- 


butions all along the chord unless shocks are present 
on both of the intermediate target pressure distribu- 
tions being averaged. If shocks are present, special 
processing (shown in fig. 7) is used for the region be- 
tween the shocks. The new shock location (*.(E)) 
will be the weighted average of the shock locations 
(x. s (B) for the forward shock and x s (C) for the aft 
shock) of the two intermediate pressure distributions. 
Thus, 

c c c 

By using the averaged pressures at the last two 
chordwise locations ahead of the forward shock, the 
new 7 pressures from the forward shock location to the 
new 7 shock location are determined by extrapolation. 
The new pressure at the foot of the shock {C p j{ E)) 
is the weighted average of pressures at the foot of 
the forward and aft shocks (C p j( B) and C p j{ C), 
respectively). Thus, 

Cpj( E)= w b c pJ (B)+(i - Wn)C pJ ( C) 

The new 7 pressures aft of the foot of the new shock 
location but upstream of the aft shock location are 
determined by interpolation. 

The new 7 target pressure distribution is used with 
the constraints to design a new 7 airfoil (either E or F) 
at the appropriate design point. (See fig. 6.) The new 
airfoil is analyzed at the other design point. These re- 
sults are used to evaluate the objective function and 
to check for any constraint violations. If the new* air- 
foil developed at one design point has a larger objec- 
tive function or constraint violations occur, the new 
target pressure distribution at the other design point 
is used to design a second new airfoil. This new 7 airfoil 
is also analyzed at the off-design point to check the 
objective function and the constraints. If the second 
new airfoil does not reduce the objective function and 
satisfy the constraints, less stringent aerodynamic 
and/or geometric constraints should be specified at 
each design point, and the target pressure-averaging 
process, including the initial processing, is repeated. 
If one of the new airfoils reduces the objective func- 
tion and satisfies the constraints, the design cycle is 
complete. To find an airfoil closer to the optimum, 
the new 7 airfoil (either E or F) becomes the baseline 
airfoil A and the target pressure- averaging procedure 
is repeated. 

Computer Resources 

For this study, the design and analysis cases for 
the dual-point design procedures were run on a Cray 
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Chart A 


Y-MP computer. One cycle of the airfoil shape aver- 
aging requires six analysis cases and two design cases 
to be run to define a new airfoil. One cycle of the tar- 
get pressure averaging requires at least five analysis 
cases and three design cases to be run. The version 
of the constrained DISC method used required only 
712 kilowords of memory for a 161 x 33 grid. Anal- 
ysis cases, which were run for typically 5500 cycles 
of the flow solver with updates for the boundary- 
layer characteristics every 100 cycles, required about 
0.028 sec of CPU time per cycle. The norm of the 
residual error was generally reduced five orders of 
magnitude from the initial free-strcam solution. De- 
sign cases were run until the maximum change in the 
chordwise pressure distribution between consecutive 
cycles was below a user-defined threshold (a change 
less than 0.004 for this study). Design cases typically 
required from 2000 to 5000 cycles of the flow solver, 
depending on the number of target pressure modifi- 
cations required by the constrained design module. 
The additional time for the computations in the con- 
strained design module led to a negligible change in 
the average CPU time required per cycle. 

Application of Revised Dual-Point 
Design Procedures 

Design Conditions 

The revised dual-point design procedures were 
used to design a new airfoil for two design points, 
one representing a climb condition and the other rep- 
resenting a cruise condition. The flow conditions for 
these points are listed in chart A. For all cases consid- 
ered, the nose-down pitching moment is constrained 
to that of the baseline airfoil so that trim drag is not a 
factor. Also, the new airfoil is constrained to have the 
same leading-edge radius and maximum thickness as 
the baseline airfoil. Different goals for the dual-point 
design procedures are selected for the two cases to 
be studied. For the first case, the selected goal was 
to reduce the combined wave drag at the two design 
points (i.e., equal weighting of the points). For the 
second case, the goal was to reduce the drag at the 
cruise condition without increasing the wave drag at 
the climb condition. 

The supercritical airfoil used as the baseline airfoil 
in reference 7 was also used as the baseline airfoil in 
this study. Design coordinates of this airfoil are listed 
in table 1. This airfoil has a maximum thickness 
of 0.115c at the 0.307c station and a leading-edge 
radius of 0.016c. The aerodynamic characteristics of 
the baseline airfoil were computed using the GAUSS2 
flow solver with the interacted boundary layer. All 
calculations were performed with the boundary- layer 


Flight condition 

Climb 

Cruise 

M x • • • • 

0.685 

0.735 

Rc 

8.9 x 10 6 

8.9 x 10 6 

Cl 

0.850 

0.608 


transition fixed on both surfaces at 0.05c. The 
u C-type” computational grid had 161 chordwise and 
33 normal grid points. The far-field boundary was 
an ellipse extending seven chords in the upstream 
direction and six chords in the normal direction from 
the trailing edge. The wake extended downstream 
6 chords from the trailing edge, and 14 chordwise 
points were in the wake. The effect of the airfoil 
lift was included in the far-field boundary conditions, 
thus allowing the far-field boundary condition to be 
brought closer to the airfoil. Computations using this 
flow solver, presented in reference 11, indicate that 
this grid is adequate for the required computations. 

The baseline airfoil was analyzed at these two 
design points, and the results are presented in fig- 
ure 8. A summary of the geometric and aerodynamic 
characteristics is found in table 2. Note that a sig- 
nificant amount of wave drag occurs at each design 
point. Two intermediate airfoils were designed, and 
the aerodynamic characteristics of each airfoil at the 
two design points are presented in figures 9 and 10. 
The wave drag of each intermediate airfoil was re- 
duced at its design point and increased at the other 
design point. Also, an increase occurred in the nose- 
down pitching-moment coefficient for both inter- 
mediate airfoils at the off-design flow condition. A 
comparison of the intermediate airfoils designed for 
the climb (airfoil B) and cruise (airfoil C) conditions 
with the baseline airfoil is presented in figure 11. The 
leading-edge radius and the maximum thickness of 
both intermediate airfoils are very close to the values 
of the baseline airfoil. Except for the forward portion 
of the airfoil, the changes in the airfoil shapes from 
the baseline for the two intermediate airfoil shapes 
are generally opposite each other. 

First Case 

In some situations, one of the design points will 
dominate the design process. Reducing the objec- 
tive function and satisfying the constraints at the 
dominant design point will lead to an airfoil that 
also satisfies the constraints at the second design 
point. The first case is an example of that situa- 
tion in that the selected goal was to reduce the com- 
bined wave drag at the two design points. The ob- 
jective function for airfoil B, designed at the climb 
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condition, is higher than that of the baseline air- 
foil A (0.0065 versus 0.0055). However, the objec- 
tive function for airfoil C, designed at the cruise con- 
dition, is lower than that of the baseline airfoil A 
(0.0040 versus 0.0055). The pitching moment of air- 
foil C at the climb condition is only slightly larger 
than that of the baseline airfoil A. If necessary, the 
dual-point design process could be repeated using a 
less negative pitching-moment constraint to achieve 
the desired value. Thus, the intermediate airfoil C, 
designed at the cruise condition, is the new airfoil. 
For this case, only the initial processing common to 
both procedures was required to define a new air- 
foil. Because the desire was to demonstrate the pro- 
cedures, the process was terminated after only the 
first cycle. To determine if additional reductions in 
the objective function are possible, the new airfoil be- 
comes the starting airfoil and the dual-point design 
process is repeated. Note that although this is an 
improved airfoil, it is not necessarily the “optimum” 
airfoil. 

Second Case 

For the second case, the objective was to reduce 
the wave drag at the cruise condition with the con- 
straint that the wave drag and pitching moment at 
the climb condition do not increase?. The objec- 
tive function and the constraint on the wave drag 
at the climb condition are plotted as a function of 
the fraction of airfoil B in figure 12. With the as- 
sumption of a linear variation of the aerodynamic 
characteristics with the change in shape, the points 
are connected with a straight line. The value of 
the objective function for airfoil B (at lip = 1) is 
greater than that of the baseline airfoil A (cruise 
cv/.M i = 0.0055 versus c dw = 0.0024), and so it is not 
acceptable. The constraint at the climb condition is 
violated (climb c dw - 0.0038 versus c dw =0.0031) 
for airfoil C (at U p = 0), and thus it is not acce pt- 
able either. Each of the revised dual-point design 
procedures will be used to produce a new airfoil. 
As noted earlier, the different characteristics are as- 
sumed to vary linearly as the percentage of the air- 
foil shape. Averaged airfoils with less than 0.26 of 
airfoil B are eliminated because of the constraint vi- 
olation (c d lc > 0.0031). The minimum value of the 
objective function achievable without violating the 
constraints at that point is selected. The weight- 
ing factor (U p) for the dual-point, design procedures 
is 0.26. 

The airfoil shape-averaging procedure (fig. 4) was 
used to define a new airfoil (D) for the second case by 
using the value of the weighting factor obtained from 
figure 12. The aerodynamic characteristics of this 


airfoil are presented in figure 13. At the climb condi- 
tion, the wave drag and pitching moment of the new 
airfoil are slightly greater than that of the baseline 
airfoil. At the cruise condition, the wave-drag coef- 
ficient has been reduced from 0.0024 to 0.0010 with 
no change in the pitching-moment coefficient. The 
shape of the airfoil redesigned by using airfoil shape 
averaging is compared with that of the baseline air- 
foil in figure 14. The new airfoil has the same max- 
imum thickness and almost the same leading-edge 
radius. Application of just 1 cycle of the airfoil shape- 
averaging procedure produced an airfoil that reduced 
the wave drag at the cruise condition but failed to 
meet the aerodynamic constraints at the climb con- 
dition. A second cycle of the airfoil shape- averaging 
procedure was used to refine the airfoil shape to sat- 
isfy the constraints. The shape of this refined air- 
foil (D‘ 2 ) is presented in figure 14. The maximum 
thickness of the refined airfoil is the same as that of 
the baseline, but the leading-edge radius is slightly 
smaller than that of the baseline airfoil. The aero- 
dynamic characteristics of the refined airfoil are pre- 
sented in figure 15. The refined airfoil satisfies all 
the aerodynamic constraints, but to satisfy the con- 
straints, the objective function (wave drag at cruise) 
has increased from 0.0010 to 0.0011. Thus, two cy- 
cles of airfoil shape averaging successfully redesigned 
the baseline airfoil. 

The values of the objective function and the wave- 
drag constraint at the climb condition for airfoil D 
have been plotted in figure 12. Two additional air- 
foils were defined: Xi using a weighting factor of 0.50, 
and X 2 using a weighting factor of 0.75. These air- 
foils were analyzed, and the values for the objective 
and constraint functions for each have also been plot- 
ted in figure 12. The values indicate that, for this 
case, the assumption of a linear variation of objective 
and constraint functions between the two intermedi- 
ate 1 airfoils, although not exact, is reasonable for the 
airfoil shape-averaging dual-point, design procedure 
for these conditions. 

The target pressure- averaging procedure? (fig. 6) 
was used to define two new airfoils for the second 
case, one for each design point. The weighting factor 
was obtained from figure 12. The weighted target 
pressure distribution for the climb condition was used 
to produce airfoil E, and the weighted target pressure 
distribution for the cruise condition was used to 
produce airfoil F. The aerodynamic characteristics of 
these two airfoils are presented in figures 16 and 17 
for both design points. 

For airfoil E, designed at the climb condition, the 
pitching-moment coefficient was unchanged, but the 
constraint on the wave-drag coefficient was violated 
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by 0.0003 at the climb condition. Both the wave- 
drag and pitching-moment coefficients were reduced 
at the cruise condition, thus reducing the objective 
function. For airfoil F, which was designed at the 
cruise condition, both the wave-drag and pitching- 
moment coefficients increased slightly at the climb 
condition, ff needed, this small increase could be 
eliminated by a second design cycle, as was done for 
airfoil shape averaging. For a practical engineering 
problem, this would probably not be considered a 
constraint violation. 

The wave-drag and pitching-moment coefficients 
were reduced at the cruise condition, thus reducing 
the objective function. A comparison of the shapes 
of these airfoils with the baseline airfoil is presented 
in figure 18. Each of the new airfoils has nearly the 
same maximum thickness and leading-edge radius as 
the baseline airfoil. The target pressure- averaging 
procedure produced two airfoils, both of which satis- 
fied the objective function and came close to satisfy- 
ing the constraints. 

The objective and constraint functions for the 
new airfoils developed using target pressure averag- 
ing (airfoils E and F) have been plotted in figure 12. 
The results indicate that an assumed linear variation 
between the intermediate airfoils is also reasonable 
for dual-point design using target pressure averaging 
for these conditions. 

The objective function and the pitching-moment 
and wave-drag constraints for the baseline airfoil, the 
intermediate airfoils, and the redesigned airfoils (D 2 , 
E, and F) are compared in figure 19. Airfoil E had 
a smaller objective function than airfoil F. However, 
airfoil E also had a minor constraint violation. A 


trade study is required to determine the better air- 
foil. The target pressure- aver aging procedure had 
more difficulty in reaching a solution that satisfied 
the aerodynamic and geometric constraints. For this 
case, airfoil shape averaging is the preferred proce- 
dure for dual-point design. 

Concluding Remarks 

Two dual-point design procedures to reduce the 
wave drag of an airfoil at two design points have 
been proposed. The first procedure to develop a re- 
designed airfoil used a weighted average of the shape 
of two intermediate airfoils redesigned at each of the 
two design points. The second procedure used a 
weighted average of two pressure distributions ob- 
tained from an intermediate airfoil redesigned at each 
of the two design points. Each procedure was used to 
design a new airfoil with reduced drag at the cruise 
condition without increasing the drag at the climb 
condition. Two cycles of the airfoil shape-averaging 
procedure successfully designed a new airfoil that re- 
duced the objective function and satisfied the con- 
straints. One cycle of the target pressure- averaging 
procedure successfully designed two new airfoils that 
reduced the objective function. However, each airfoil 
had a minor constraint violation that could be elim- 
inated with additional design cycles. For this study, 
airfoil shape averaging was the preferred method for 
dual-point design. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
June 15, 1994 
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Appendix 

Preliminary Dual-Point Design 
Procedures 

Preliminary dual-point design procedures were 
defined in reference 7. These earlier procedures were 
applied to redesign an airfoil to reduce the drag at 
two design points. The baseline airfoil and the re- 
designed airfoil were tested in a wind tunnel to verify 
the drag reduction of the redesigned airfoil relative to 
the baseline. These results indicated that averaging 
is a viable technique for a dual-point design proce- 
dure and that the computer code used to model the 
flow field does a good job in simulating the experi- 
mental flow field. For completeness, the preliminary 
dual-point design procedures, their application to re- 
duce the drag at two design points, and results from 
the experimental verification of the improvements are 
presented in this appendix. These procedures used 
an early version of the airfoil constrained design tech- 
nique (the constrained DISC) described in the main 
section of this report, and thus a description will not 
be repeated in the appendix. 

Description of Preliminary Dual-Point 

Design Procedures 

Two preliminary dual-point design procedures 
were defined in reference 7. Each procedure used 
an iterative approach to redesign an airfoil. The first 
procedure to develop a new airfoil used a weighted 
average of the shapes of two intermediate airfoils re- 
designed at each of the two design points. This pro- 
cedure will be referred to as “airfoil shape averag- 
ing/' The second procedure used a weighted average 
of two pressure distributions obtained from an inter- 
mediate airfoil redesigned at the same design point. 
This procedure will be referred to as “target pressure 
averaging." A description of each of these dual-point 
design procedures follows. 

Airfoil shape averaging. The first dual-point 
design procedure averages the shapes of two re- 
designed intermediate airfoils to develop a new air- 
foil. A schematic of one design cycle of this proce- 
dure is presented in figure Al. The starting airfoil is 
first analyzed to obtain the chordwise pressure dis- 
tribution at each design point. These pressure dis- 
tributions are adjusted by the user to define the ini- 
tial target pressure distributions at each design point. 
The user-specified adjustments allow a specific char- 
acteristic. such as reduced shock strength or reduced 
adverse pressure gradient, to be included in the de- 
sign process. Starting from the appropriate initial 
target pressure distribution, the constrained design 
technique develops an intermediate airfoil at each de- 
sign point. A new airfoil is obtained from the simple 


average of the shapes of the two intermediate airfoils. 
The new airfoil is analyzed at both design points, and 
if the results are not satisfactory, the new airfoil be- 
comes the starting airfoil and the process is repeated. 
This procedure can be expanded to more than two 
design points by simply adding additional process- 
ing paths for additional design points and defining a 
suitable weighting function for the averaging of the 
intermediate airfoil shapes. 

Airfoil target pressure averaging. To de- 
velop an improved airfoil, the second dual-point de- 
sign procedure averages the detailed target pressure 
distributions developed in the iterative process rather 
than the airfoil shapes. A schematic of the initializa- 
tion and one design cycle is presented in figure A2. 
The starting airfoil is first analyzed to obtain the 
chordwise pressure distributions at each of the design 
points. The pressure distributions from the analysis 
are adjusted to define the initial target pressure dis- 
tributions at each design point. The initialization 
can start with either design point, and subsequent 
design cycles will use alternating design points. For 
this example, design point 1 was selected to start 
the process. The adjusted target pressure distribu- 
tion at design point 1 is used with the constraints to 
design the first intermediate airfoil. Tin' analysis of 
the starting airfoil, the definition of the init ial target 
pressure distributions, and the design of the first in- 
termediate airfoil constitute the initialization of the 
process. 

The initialization is followed by multiple airfoil 
design cycles. The latest intermediate airfoil is an- 
alyzed to obtain the pressure distribution at each 
design point. These pressure distributions are used 
as intermediate target pressure distributions for each 
design point. The last two intermediate target pres- 
sure distributions for the design point are averaged 
to define the next target pressure distribution, and 
this averaged distribution is t hen used with the con- 
straints to design the next intermediate airfoil. The 
analysis of an intermediate airfoil, the averaging of 
the last two intermediate target pressure distribu- 
tions, and the design of the next intermediate airfoil 
constitute one airfoil design cycle. The intermedi- 
ate airfoil is analyzed at both design points. If the 
results are not satisfactory, the airfoil design cycle 
is repeated at the other design point to develop the 
next intermediate airfoil. The airfoil design cycle is 
repeated at alternate design points as needed. 

Application of Preliminary Dual-Point 

Design Procedures 

A supercritical airfoil shape was selected as 
the baseline airfoil (A) for this study, and design 


10 



coordinates are listed in table 1. This airfoil has a 
maximum thickness of 0.115c at the 0.307c station 
and a leading-edge radius of 0.016c. The baseline 
airfoil aerodynamic characteristics were computed 
using the GAUSS2 flow solver wuth the interacted 
boundary layer. All calculations were performed 
with the boundary- layer transition fixed on both sur- 
faces at 0.05c. Two design points were selected: 
one representing a long-range cruise condition and 
the other representing a high-speed cruise condition. 
The flow conditions for these points are listed in 
chart B. Results from the analysis of the baseline air- 
foil at these design points are presented in figure A3, 
and a summary of the geometric and aerodynamic 
characteristics is found in table 3. 


Chart B 


Cruise condition 

Long range 

High speed 

Moc • ■ • • 

0.654 

0.735 

Rr 

4.5 x 10 r> 

8.9 x 10 6 

( 'l 

0.979 

0.508 


The preliminary versions of the procedures for 
airfoil shape averaging and target pressure averag- 
ing were applied in an attempt to redesign the base- 
line airfoil. The objective function was the sum 
of the wave drag at the two design points. The 
section lift, wave drag, and pitching-moment coeffi- 
cients from the baseline airfoil were used as the aero- 
dynamic constraints on the design of the improved 
airfoils. The pitching-moment-coefficient constraint 
was a one-sided constraint in that the pitching- 
moment coefficient of the new airfoil could not be 
more negative than the pitching-moment coefficient 
of the baseline airfoil. Similarly, the wave-drag co- 
efficient of the new airfoil could not be greater than 
the wave drag of the baseline airfoil. The leading- 
edge radius, maximum thickness, and thickness at 
the 0.85c location of the baseline airfoil were used as 
the geometric constraints. In general, the constraint 
on the local thickness at 0.85c did not impact the de- 
sign process, and thus it will not be discussed herein. 

By using airfoil shape averaging, an intermediate 
airfoil was designed at each design point and was 
analyzed at the other design point. The aerodynamic 
results (presented in figs. A4 and A5) indicate that 
even though each intermediate airfoil had less wave 
drag than the baseline airfoil at the particular design 
point, it had a higher wave drag than the baseline 
airfoil at the other design point, thus violating one of 
the constraints. The value of the objective function 
for each of the intermediate airfoils was larger than 


that of the baseline airfoil. A comparison of the two 
intermediate airfoil shapes with the baseline airfoil 
shape (presented in fig. A6) shows that changes in 
shape for the two intermediate airfoils relative to the 
baseline were generally in opposition to each other. 

The preliminary procedures of airfoil shape aver- 
aging and target pressure averaging that used simple 
averaging could not design a new airfoil that satisfied 
all the constraints and that significantly reduced the 
sum of the wave-drag coefficients at the two design 
points. To develop a new airfoil, the design space 
was temporarily expanded to see if an alternate path 
could be found around the boundary posed by the 
constraints. The preliminary version of the target 
pressure-averaging procedure was used with relaxed 
constraints to develop an airfoil with reduced wave 
drag at both design points. The preliminary ver- 
sion of the airfoil shape-averaging procedure was then 
used to refine the shape to satisfy all constraints. 

Target pressure averaging was used wfith relaxed 
constraints to design a series of intermediate air- 
foils. The leading-edge-radius constraint was elim- 
inated and the pitching-moment constraint was re- 
laxed to -0.077 for the long-range cruise point. The 
design value for the wave-drag coefficient w r as 0.0010 
for the long-range cruise point and 0.0001 for the 
high-speed cruise point. The baseline pressure dis- 
tributions from figure A3 were modified to reduce 
the Mach number (negative pressure coefficient) just 
ahead of the shock. The modified pressure distribu- 
tions became the initial target pressure distributions 
for the target pressure-averaging procedure. Sev- 
eral iterations of the preliminary target pressure- 
averaging procedure were used to define a series of 
intermediate airfoils. The fifth intermediate airfoil 
showed a significant reduction in the wave drag at the 
long-range cruise design point without a significant 
increase at the high-speed cruise design point. This 
airfoil satisfied the maximum thickness constraint. 
However, the design process reduced the leading- 
edge radius from 0.016c to 0.005c. Also, the negative' 
pitching moment exceeded the constrained value at 
both design points. This airfoil must be redesigned 
to increase the leading-edge radius and to reduce the 
negative pi telling- moment coefficient. 

The airfoil shape- averaging procedure was used 
with the full set of constraints to redesign the air- 
foil with the reduced leading-edge radius. The aero- 
dynamic characteristics of the redesigned airfoil are 
presented in figure A7. Relative to the baseline 
airfoil, wave drag has been reduced at both design 
points. The pitching moment was close to the con- 
strained value. At the long-range cruise point, the 
Mach number ahead of the shock has been reduced 
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and the loading on the forward portion of the lower 
surface has been increased. 

Figure A8 presents a comparison of the shapes of 
the baseline airfoil arid redesigned airfoil. The max- 
imum thickness and leading-edge radius of the re- 
designed airfoil are close to the values of the baseline 
airfoil. Coordinates of the redesigned airfoil are pre- 
sented in table 4. Through the sequential application 
of the two preliminary dual-point design procedures, 
a redesigned airfoil was developed that, for practical 
purposes, does not have any violations of the geo- 
metric or aerodynamic constraints and that reduces 
the wave drag at both design points. 

Experimental Verification of Redesigned 

Airfoil 

Two airfoil models were built, one incorporating 
the baseline airfoil shape and the other incorporating 
the redesigned airfoil shape. Each model was tested 
in the Langley 0. 3-metcr TCT with the 13- by 13-in. 
two-dimensional adaptive-wall test section installed 
in the circuit. (See ref. 8.) The 0.3-m TCT is a fan- 
driven, cryogenic pressure tunnel that uses nitrogen 
as a test gas. It can test airfoil models over a 
Mach number range from about 0.20 to 0.95 at chord 
Reynolds numbers up to about 50 x 10°, based on a 
model chord of 6 in. 

The models used in these tests had a 6-in. chord, 
and both airfoil models were tested with transition 
strips placed at the 5-percent-chord location on both 
surfaces of the model. Details of the experimental 
tests of the baseline airfoil model are reported in ref- 
erence 9. A comparison of the experimental and cal- 
culated chord wise pressure distributions is presented 
in figure A9 for the baseline airfoil and in figure A 10 
for the redesigned airfoil. In general, the pressures 
are in good agreement. However, the calculated 
shock location is usually upstream of the measured 
location and the calculated trailing-edge pressure is 
more positive than the measured value. This can be 
attributed to the very simple boundary-layer model 


used for this study. Overall, the agreement indicates 
that the flow solver with the simple boundary layer 
does a good job in simulating the experimental flow 
field. 

A comparison of the measured chordwise pressure 
distributions for both airfoils at test conditions close 
to the two design points is presented in figure All. 
At both design points, the loading has been increased 
on the lower surface over the front portion of the re- 
designed airfoil, and the negative pressure coefficient 
on the upper surface just ahead of the shock has been 
reduced. These changes are similar to those found 
in the computed results presented in figure A 7 . A 
comparison of the experimental integrated force and 
moment coefficients for the two airfoils is presented 
in figure A 12. The redesigned airfoil reduced the 
drag coefficient by 0.0022 at the long-range cruise 
condition (A/^ = 0.654) and reduced the drag co- 
efficient by 0.0005 at the high speed cruise condition 
(Moc — 0.735). At each design point, little difference 
occurs in the pitching-moment coefficients between 
the baseline and the redesigned airfoils. Because the 
drag reductions at A/ 0 c = 0.654 are even larger at val- 
ues of (■[ above the design value, the redesigned airfoil 
probably has a larger buffet margin than the base- 
line airfoil. These drag reductions include not only 
changes in the wave drag but also changes in the 
skin friction and pressure drag. At both the long- 
range cruise point and the high-speed cruise point, 
the experimental reduction in the total drag co- 
efficient (fig. A 12) was slightly larger than the 
calculated reduction in the wave-drag coefficient 
(fig. A7). 

The agreement of the experimental results and 
the calculations for the pressure coefficients indicates 
that the flow solver with the interacted boundary 
layer is suitable for modeling the flow fields for the 
design conditions studied. The agreement of the pre- 
dicted wave-drag reduction and the measured total 
drag reduction indicates that averaging is a viable 
concept for a dual-point design procedure. 
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Figure Al. Preliminary version of multipoint design procedure using airfoil shape averaging. 
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Figure A 12. Comparison of experimental force and moment coefficients of baseline and redesigned airfoils 
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Table 1. Design Coordinates of Baseline Airfoil 
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Upper surface 

Lower surface 


x/c 

Upper surface 

Lower surface 

0.00000 

0.00000 

0.00000 


0.44557 

0.06691 

-0.04100 

0.00099 

0.00635 

-0.00489 


0.46597 

0.06601 

-0.03958 

0.00301 

0.01117 

-0.00821 


0.48646 

0.06488 

-0.03808 

0.00604 

0.01562 

-0.01132 


0.50699 

0.06353 

-0.03651 

0.01005 

0.01974 

-0.01431 


0.52756 

0.06197 

-0.03487 

0.01500 

0.02362 

-0.01702 


0.54812 

0.06020 

-0.03317 

0.02088 

0.02731 

-0.01949 


0.56865 

0.05826 

-0.03141 

0.02764 

0.03076 

-0.02183 


0.58912 

0.05617 

-0.02961 

0.03528 

0.03395 

-0.02407 


0.60950 

0.05397 

-0.02777 

0.04374 

0.03692 

-0.02622 


0.62977 

0.05168 

-0.02591 

0.05302 

0.03969 

-0.02830 


0.64990 

0.04933 

-0.02403 

0.06308 

0.04230 

-0.03035 


0.66986 

0.04692 

-0.02214 

0.07389 

0.04477 

-0.03234 


0.68962 

0.04448 

-0.02027 

0.08543 

0.04713 

-0.03428 


0.70915 

0.04200 

-0.01842 

0.09766 

0.04937 

-0.03617 


0.72843 

0.03948 

-0.01662 

0.11056 

0.05152 

-0.03797 


0.74742 

0.03694 

-0.01489 

0.12411 

0.05358 

-0.03968 


0.76611 

0.03438 

-0.01324 

0.13826 

0.05554 

-0.04126 


0.78445 

0.03181 

-0.01170 

0.15300 

0.05740 

-0.04270 


0.80243 

0.02922 

-0.01028 

0.16830 

0.05915 

-0.04400 


0.82002 

0.02665 

-0.00897 

0.18413 

0.06078 

-0.04512 


0.83718 

0.02409 

-0.00781 

0.20045 

0.06228 

-0.04605 


0.85389 

0.02157 

-0.00678 

0.21725 

0.06364 

-0.04680 


0.87013 

0.01910 

-0.00591 

0.23450 

0.06484 

-0.04735 


0.88585 

0.01670 

-0.00520 

0.25216 

0.06587 

-0.04769 


0.90105 

0.01438 

-0.00463 

0.27021 

0.06674 

-0.04783 


0.91568 

0.01217 

-0.00423 

0.28863 

0.06743 

-0.04777 


0.92972 

0.01006 

-0.00397 

0.30737 

0.06796 

-0.04751 


0.94314 

0.00807 

-0.00385 

0.32642 

0.06831 

-0.04705 


0.95592 

0.00621 

-0.00386 

0.34575 

0.06851 

-0.04642 


0.96802 

0.00447 

-0.00398 

0.36533 

0.06854 

-0.04561 


0.97942 

0.00285 

-0.00421 

0.38513 

0.06840 

-0.04465 


0.99009 

0.00136 

-0.00453 

0.40512 

0.06809 

-0.04355 


1.00000 

0.00000 

-0.00490 

0.42527 

0.06760 

-0.04233 





28 





Table 2. Geometric and Aerodynamic Characteristics of 
Baseline Airfoil at Climb and Cruise Design Conditions 


Design conditions 

Climb 

Cruise 

A/qc .... 

0.685 

0.735 

R c 

8.9 x 10 6 

8.9 x 10 6 

q, deg ... 

2.88 

1.17 

Cl 

0.850 

0.608 

c d,w * ■ * ■ 

0.0031 

0.0024 

Cm 

-0.060 

-0.079 

x s /c .... 

0.35 

0.48 

rijc .... 

0.016 

0.016 

tmax/c • • • 

0.115 

0.115 


Table 3. Geometric and Aerodynamic Characteristics of Baseline 
Airfoil at Long-Range and High-Speed Cruise Design Conditions 


Cruise conditions 

Long range 

High speed 

Mqo .... 

0.654 

0.735 1 

R c 

4.5 x 10 6 

8.9 x 10 6 

q, deg ... 

4.37 

0.69 

q 

0.979 

0.508 


0.0042 

0.0007 

Cm 

-0.054 

-0.077 

x s /c .... 

0.23 

0.44 

r le /c .... 

0.016 

0.016 

^max/c • ■ • 

0.115 

0.115 




Table 4. Design Coordinates of Redesigned Airfoil 



x/c 

Values of y/c for 


x/c 

Values of y/c for 

Upper surface 

Lower surface 


Upper surface 

Lower surface 

0.00000 

0.00000 

0.00000 


0.44557 

0.06430 

-0.04819 

0.00099 

0.00619 

-0.00601 


0.46597 

0.06390 

-0.04692 

0.00301 

0.01045 

-0.00987 


0.48646 

0.06331 

-0.04541 

0.00604 

0.01431 

-0.01315 


0.50699 

0.06251 

-0.04369 

0.01005 

0.01788 

-0.01594 


0.52756 

0.06147 

-0.04178 

0.01500 

0.02122 

-0.01832 


0.54812 

0.06018 

-0.03969 

0.02088 

0.02438 

-0.02035 


0.56865 

0.05865 

-0.03744 

0.02764 

0.02740 

-0.02210 


0.58912 

0.05689 

-0.03506 

0.03528 

0.03030 

-0.02362 


0.60950 

0.05492 

-0.03257 

0.04374 

0.03307 

-0.02496 


0.62977 

0.05277 

-0.02999 

0.05302 

0.03570 

-0.02617 


0.64990 

0.05046 

-0.02735 

0.06308 

0.03817 

-0.02732 


0.66986 

0.04802 

-0.02469 

0.07389 

0.04048 

-0.02847 


0.68962 

0.04546 

-0.02204 

0.08543 

0.04269 

-0.02967 


0.70915 

0.04282 

-0.01945 

0.09766 

0.04480 

-0.03096 


0.72843 

0.04012 

-0.01695 

0.11056 

0.04685 

-0.03235 


0.74742 

0.03738 

-0.01459 

0.12411 

0.04881 

-0.03386 


0.76611 

0.03462 

-0.01240 

0.13826 

0.05068 

-0.03546 


0.78445 

0.03186 

-0.01042 

0.15300 

0.05248 

-0.03714 


0.80243 

0.02913 

-0.00867 

0.16830 

0.05419 

-0.03888 


0.82002 

0.02643 

-0.00718 

0.18413 

0.05585 

-0.04063 


0.83718 

0.02380 

-0.00594 

0.20045 

0.05743 

-0.04235 


0.85389 

0.02125 

-0.00496 

0.21725 

0.05891 

-0.04400 


0.87013 

0.01880 

-0.00423 

0.23450 

0.06025 

-0.04554 


0.88585 

0.01646 

-0.00372 

0.25216 

0.06145 

-0.04692 


0.90105 

0.01423 

-0.00342 

0.27021 

0.06246 

-0.04812 


0.91568 

0.01212 

-0.00329 

0.28863 

0.06328 

-0.04911 


0.92972 

0.01011 

-0.00330 

0.30737 

0.06391 

-0.04988 


0.94314 

0.00822 

-0.00343 

0.32642 

0.06437 

-0.05041 


0.95592 

0.00641 

-0.00365 

0.34575 

0.06466 

-0.05069 


0.96802 

0.00469 

-0.00392 

0.36533 

0.06481 

-0.05072 


0.97942 

0.00305 

-0.00423 

0.38513 

0.06483 

-0.05049 


0.99009 

0.00149 

-0.00456 

0.40512 

0.06475 

-0.04999 


1.00000 

0.00000 

-0.00489 

0.42o2 7 

0.06458 

-0.04922 
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Initial target 
pressure distribution 



Figure 1. Processing for constrained-design technique of airfoil 











Original airfoil surface 



xx 


(a) Modification of surface shape for curvature change. 



Modified airfoil surface -* 

(b) Surface shearing to close trailing edge. 

Figure 3. Modification of airfoil shape for change in surface curvature 



Initial processing I Airfoil shape averaging: 

I D - W B B + (1 - W B )C 
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Figure 4. Multipoint design procedure using airfoil shape averaging. 



Objective function: F 

ith constraint: G < G max Constraint 

violation: 
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Figure 5. Determination of weighting factor (W) for objective function (F) with constraint function (G). 
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Figure 6. Multipoint design procedure using target pressure averaging. 




_ Upstream of 

Upstream of 


(a) 

both shocks 

averaged shock 



Weighted average 

Extrapolation 

c pj 

(b) 


Downstream of 
both shocks 


Weighted average 


Location 


Cp computation 


(a) : Downstream of averaged shock 

(b) : Interpolation 

Figure 7. Weighted averaging of target pressure distributions when shocks are present. 
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Figure 12. Determination of weighting factor for multipoint procedures. Solid symbols with letter denote new 
airfoils. 
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(b) Cruise. = 0.735; q - 0.608. 


Figure 16. Comparison of analysis of baseline airfoil A with analysis of airfoil E (developed using target c 
averaging at climb design point). 
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Figure 19. Comparison of objective function and constraint functions of baseline, intermediate, and now airfoils 
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